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ABSTRACT 

We revisit and extend the analysis supporting a 60 year-old suggestion that cosmic 
rays air showers resulting from primary particles with energies above 10^^ eV should 
be straightforward to detect with radar ranging techniques, where the radar echoes are 
produced by scattering from the column of ionized air produced by the shower. The 
idea has remained curiously untested since it was proposed, but if our analysis is correct, 
such techniques could provide a significant alternative approach to air shower detection 
in a standalone array with high duty cycle, and might provide highly complementary 
measurements of air showers detected in existing and planned ground arrays such as the 
Fly's Eye or the Auger Project. The method should be particularly sensitive to showers 
that are transverse to and relatively distant from the detector, and is thus effective in 
characterizing penetrating horizontal showers such as those that might be induced by 
ultra-high energy neutrino primaries. 



Subject headings: Cosmic-rays, Cosmic-ray detectors. Extensive air showers, Neutrinos, 
Radar 
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1. Introduction 

Extensive air showers (EAS) resulting from 
primary cosmic-ray particles of energies above 1 
EeV (10^^ cV) produce an ionization trail which 
is comparable to that of micro-meteors, which 
have been detected for many decades using radar 
methods (Lovell 1948; Greenhow, 1952; Hanbury 
Brown k Lovcll 1962). In fact, Blackctt & Lovcll 
(1940) proposed that earlier detections of spo- 
radic radio reflections at altitudes of 5-20 km 
from the troposphere and mesosphere (cf. Col- 
well & Friend 1936; Watson Watt et al. 1936a, 
1936b, 1937; Appleton & Piddington 1938) might 
be attributed to cosmic ray EAS, although at the 
time both the energies and fluxes of such events 
were only crudely known. Later confirmation of 
the meteoric origin of higher altitude ionization 
trails may have contributed to a lack of follow- 
up on this prescient suggestion for the source of 
these low-altitude sporadics. The idea was re- 
examined in the early 1960's (Suga 1962) and a 
proposed experiment was described by Matano 
et al. (1968); however, no results from this test 
have appeared in the literature to date. 

EAS ionization trails are now commonly de- 
tected by their air fluorescence emission at near- 
UV wavelengths (cf. Baltrusaitas et al. 1985). 
Future large EAS detector arrays such as the 
Auger project (Guerard et al. 1998) and the pro- 
posed space mission OWL/AirWatch^ (Scarsi et 
al. 1999; Krizmanic et al. 1999) have made fluo- 
rescence detection of EAS a centerpiece of their 
approach, since it can provide information such 
as the position of the shower maximum and the 
total shower energy which can be difficult to pin 
down with particle detectors alone. 

One of the most compelling reasons to ex- 
tend the sensitivity of EAS detectors in this en- 
ergy regime is the possibility that neutrinos of 

^Recently, the Air Watch program has evolved into a mis- 
sion under consideration for the International Space Sta- 
tion, and has been renamed Extreme Universe Space Ob- 
servatory (EUSO). 



energy > 10 eV may be an important com- 
ponent of the primary particles (Capelle et al. 
1999), contributing to the fluxes at these ener- 
gies either directly or indirectly. Neutrinos and 
other highly penetrating particles will induce air 
showers that develop much deeper in the atmo- 
sphere than hadron-initiated showers. A unique 
signature of their presence would be the detec- 
tion of highly inclined (near-horizontal) showers 
for which the point of maximum development is 
many radiation lengths deeper than is possible 
for a hadronic or electromagnetic shower. Thus 
any technique which (like the air fluorescence ap- 
proach) has a particular sensitivity to showers of 
this type is of interest in addressing this problem. 

To our knowledge no one has reported any 
studies the ionization trail of EAS using radar 
echo techniques, although as we will show here, 
the signals should be clearly detectable using 
standard radar methods. One reason for this may 
be the transient nature of the expected echoes; 
short duration transients are typically rejected by 
radar systems designed to detect or track targets 
where the reflectivity does not rapidly decay. 

Ionization columns: meteors vs. EAS. 

Although there are many significant differences 
between the ionization caused by meteors and 
that of EAS, we can use the basic formalism de- 
veloped for meteor study as a starting-point for 
analysis of EAS ionization. Meteor ionization 
trails are commonly parameterized in terms of 
their ionization line density a (electrons m~^), 
a measure of the total ionization content divided 
by the length of the meteor track. Typical radar- 
detected meteors occur at heights of 80-120 km, 
and have line densities of a = 10^^ to 10^^ m~^. 
At the lowest detectable line densities, the inci- 
dent meteor has a mass of ~ 1 fig, with a radius 
less than 100 /im. At typical velocities of > 30 
km s^^, the implied kinetic energy of these me- 
teor grains is 0.05 Joules or more, much of which 
goes into ionization of the air along it path. 
A cosmic ray proton of energy 10^^ eV also 



2 



has a kinetic energy of order 0.1 J, and much of 
this energy also ultimately ends up in the form of 
ionization and excitation of atoms of the air along 
the path of the shower of charged particles that 
results from the proton's collision with nuclear 
hadrons. The primary differences between the 
meteor track and cosmic-ray-induced EAS are in 
the way the ionization column forms, and in the 
resulting ionization density profile. 

For the meteor, ablation of material from 
its surface yields atoms with kinetic energies of 
10^ - 10^ eV, which ionize air molecules by di- 
rect collision with a mean free path of several 
cm. This produces an ionization column with 
an approximately uniform distribution of radial 
density, and an initial radius of 2-10 m. The 
density then evolves with time due to diffusion, 
convective processes, bulk motions of the air, and 
the Lorentz forces of the ambient electric and ge- 
omagnetic fields. Electron attachment and re- 
combination eventually complete the process of 
dissipating the ionization column. At radar fre- 
quencies in the HF to lower VHF range (10-100 
MHz), echoes from typical meteors may be de- 
tectable for several seconds after the meteor is 
gone. 

The ionization in an EAS, in contrast, is not 
produced by a single body, but rather by the col- 
lective effects of the disk highly energetic parti- 
cles (mostly electrons and positrons) that make 
up the body of the shower. Because the lat- 
eral distribution of these particles spreads out as 
the shower progresses, the ionization column is 
formed with a different initial distribution than 
that of a meteor, reflecting the evolution of the 
cross-sectional charged particle density. In addi- 
tion, since the shower propagates essentially at 
the speed of light, it appears almost instanta- 
neously compared to even the fastest meteors at 
~ 100 km s-^ 

An important measure of the transverse charged 
particle distribution in an EAS is the Moliere ra- 
dius Vm within which of order 90% of the charged 
particles can be found. For air at sea level. 



Tm — 70 m, but it is important to note that 
within r„i the radial distribution is a power law, 
and most showers retain a tight core of parti- 
cles of diameter several m or less which may look 
much like the column initially produced by a me- 
teor. In fact, of order 10% of the total number of 
electrons in an EAS are likely to be found within 
a radius of ~ rmf^O, comprising much less than 
1% of the area of the Moliere disk. 

Lightning ionization columns. Lightning dis- 
charges also produce intensely ionized regions 
which have also been the subject of many radar 
studies since the early 1950's (Ligda 1950, 1956; 
Miles 1952). Ionization charge densities pro- 
duced by lightning are many orders of magnitude 
higher than those produced by either EAS or 
small meteors, and the plasma channels produced 
by lightning are also of much smaller diameter 
than either meteor or EAS ionization columns. 
However they do occasionally occur in clustered 
networks with transverse scales comparable to 
those of EAS and meteors, and at altitudes that 
overlap the altitudes of EAS maxima. Thus they 
can provide examples of radar targets with fea- 
tures that are in some ways relevant to our dis- 
cussion. 

In the following section we present the con- 
cept of the radar cross section and outline specific 
cases relevant to EAS ionization columns. Sec- 
tion 3 develops a semi-analytical method for de- 
termining the ionization densities for a given air 
shower energy and altitude, and discusses some 
issues related to the evolution of the electron den- 
sity. Section 4 then makes estimates of the radar 
echo power based on the formalism developed in 
the previous two sections. Section 5 concludes 
the paper with a discussion of some applications 
to specific experimental conditions. 
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2. Radar detection of ionization columns 

For radar detection of the columns that result 
from meteors, EAS events, or lightning, there 
are two regimes to consider, depending on the 
plasma frequency Vp of the ionized region: 

Up = ^ {uec'^ /nme) 8.98 x 10^ Hz (1) 

where ng is the electron density in cm~^. These 
two regimes are known in as the under- and over- 
dense regimes, respectively, and traditionally (in 
the meteor radar literature) arc divided at the 
line density of a ~ 10^^ m~^. However, because 
of the significant differences in the radial distri- 
bution of the electrons in EAS ionization com- 
pared to that of meteors, we distinguish them 
only on the basis of the ratio of radar frequency 
to plasma frequency, a distinction that is also 
common to radar lightning measurements. Thus 
we will consider the over-dense portion of an 
ionization column to be that region where the 
electron density is high enough that the plasma 
frequency exceeds the radar frequency, and the 
radar cannot penetrate it and is reflected. Un- 
derdense columns are those which have electron 
densities such that the local plasma frequency is 
below the frequency of the incoming radar, which 
can therefore penetrate the ionized region. ^ 

Radar targets are most commonly described 
in terms of their effective radar cross section 
(RCS), a measure of the equivalent physical area 
of an ideal scattering surface. Here we will de- 
velop concepts useful to understanding the RCS 
characteristic to ionization columns, to prepare 
for later estimation of the expected radar return 
power from EAS-induced ionization. 



^We stress here that underdense ionization columns may 
still produce a large radar return, for although the incident 
radio waves penetrate the column, the scattering may still 
be largely coherent, yielding a return power that scales 
quadratically with the electron density. This feature is 
of particular importance in determining radar's ability to 
provide estimates of the shower electron density and thus 
ultimately the energy of the primary particle. 



In the following subsections we introduce gen- 
eral formulas for the RCS for both over- and 
under-dense cases. An additional distinction is 
also useful in each case: the long wavelength, or 
Rayleigh regime involves targets with characteris- 
tic dimensions smaller than the radar wavelength 
used; and the short wavelength, or optical regime 
involves cases where the characteristic target size 
is much larger than the radar wavelength. 

2.1. Overdense ionization columns 

When the electron density ne(rc) is high enough 
to produce a surface where the plasma frequency 
exceeds the frequency of the incident radiation at 
some critical radius from the track, the result- 
ing index of refraction becomes imaginary, and 
total external reflection of the radiation occurs. 
Under these conditions the surface at Vc can be 
treated to flrst order as a metal cylinder, and the 
RCS is accordingly greatly enhanced. 

2.1.1. Overdense columns: Rayleigh regime 

When Tc ^ A, the cylindrical region of re- 
flecting plasma can be approximated as a thin 
wire. As wc will show later when we develop 
an estimate of the ionization density for EAS, 
this is the case that will be most relevant to low- 
frequency radar echoes from EAS. In this case 
the maximum of the RCS, which occurs for nor- 
mal incidence, can be approximated as (Crispin 
& Maffett 1965): 
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where (f> the angle of linear polarization with re- 
spect to the axis of the wire. The behavior of 
ab at angles other than normal incidence is given 
by: 
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where is the angle as measured from the axis 
of the wire (or EAS in our case). This approxi- 
mation is vahd in the range 60° < 9 < 120°, but 
underestimates the cross section at angles closer 
to the axis of the wire, due to what arc known 
as traveling wave effects, which produce large 
cross sections near the axis, in some cases com- 
parable to that at normal incidence (Medgyeshi- 
Mitscliang & Putnam 1985), and significantly 
broader in angular range. However, the treat- 
ment of such effects is very complex and beyond 
our scope at present. 

2.1.2. Overdense columns: Optical regime 

When rc 3> A, the reflection from the critical 
surface becomes almost entirely specular, and for 
normal incidence we have (Kraus 1988) 



27rrc-L2 
A 



(4) 



where ah^max denotes the maximum value of the 
RCS, and L is the length of the cylindrical reflect- 
ing surface. Radar cross sections for metal cylin- 
ders have been studied in detail (cf. Mcdgycshi- 
Mitschang &; Putnam 1985) and are complicated 
functions of angle, with resonances and nulls that 
depend strongly on the length and other fac- 
tors. As we shall establish below, EAS ioniza- 
tion columns will generally not fall into the opti- 
cal regime for the overdense case, except at very 
high EAS energies and very low radar frequen- 
cies. For this reason we do not treat these effects 
further here. 

2.2. Underdense ionization columns 

For radar frequencies well above the plasma 

frequency Up, the radiation penetrates the entire 
ionization column, and the electrons within the 
ionization column scatter independently accord- 
ing to the Thomson cross section 



Svr 



mpC^ 



= 6.65 X 10-^^m^ . (5) 



The total effective radar backscatter cross sec- 
tion CTfe will then depend on the individual phase 
factors of each of the scattering electrons. 

2.2.1. Underdense columns: Rayleigh regime 

For Tm A, and angles that are nearly per- 
pendicular to the track, all of the electrons scat- 
ter coherently over a longitudinal region of the 
track Lp = yJ\R/2 where R is the perpendicular 
distance to the track and A the radar wavelength 
(Hanbury-Brown & Lovell 1962). This region of 
the track is known as the first Fresnel zone of 
the track, and the total radar backscatter cross 
section then becomes 



(6) 



where = aLp = a^J XR/2 is the number of 
electrons within a single Fresnel zone along the 
track. The radar cross section depends of the 
square of the electron density because of the as- 
sumption of full coherence. 

Except for radar frequencies / <C 10 MHz, 
most radar ranging of EAS ionization columns is 
not in the Rayleigh regime. However, the concept 
of the Fresnel zone, which establishes a charac- 
teristic length scale over the longitudinal extent 
of the track, will still be useful in what follows. 

2.2.2. Underdense columns: Optical regime 

For tracks where the radar frequency f > i^p, 
but 

^ A/4, the assumption of coherent scat- 
tering is no longer satisfied, and the individual 
phase factors of the electrons must be included. 
For this case the effective RCS can be written 
(cf. Wehner 1987; also see appendix): 



(7) 



where k is the wave vector of the incident field 
(k = 2tt/X), r is the vector distance to the vol- 
ume element at which the scattering takes place. 
Here we are neglecting refractive effects of the 
ionization column on the incident radiation, and 
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assuming that the incident and backscattered 

waves satisfy k • r ~ \kr\. 

Note that the argument of the exponential 
here includes an extra factor of two to represent 
the two-way phase. This is due to the fact that as 
the scattered radiation returns on the same path 
as the incident radiation, it picks up an addi- 
tional phase factor equal to that of the incoming 
radiation. The term ne(r)y^ or d^r represents 
the differential contribution to the scattered elec- 
tric field of a volume element of electrons which 
scatter coherently. 

If we now write q = 2k then equation 7 be- 
comes 



^^(q) 



(8) 



Equation 8 thus reduces the problem of es- 
timating the effective cross section in the un- 
dcrdensc case to that of calculating the Fourier 
transform (and the resulting power spectrum) of 
the electron number density distribution. In the 
appendix we describe the geometry more explic- 
itly and show how this relation arises. 

A broader implication of this result can be 
stated as follows: A measurement of the com- 
plex amplitude of the radar echo from an EAS 
in the underdense case is proportional to a mea- 
surement of one Fourier component of the shower 
ionization profile. This will be true even in the 
case of a bistatic radar (receiver location not co- 
incident with transmitter) although there will 
be an additional coefficient to account for the 
bistatic cross section and its angular behavior. 

3. EAS ionization densities 

In the previous section we introduced the 
mechanisms for radar reflection from cylindrical 
plasma columns in two regimes of ionization den- 
sity. Here we estimate the spatial distribution 
of the expected ionization from a high energy 
cosmic-ray air shower, which will then allow us 
to assess their effective radar cross section. 



3.1. Longitudinal ionization 

There are many years of development of the 

theory of EAS production. The most accurate 
treatments of the evolution of the electron den- 
sity in the shower require numerical simulations, 
but there are a number of parameterizations 
available that yield results accurate enough for 
our needs. Here we use the analytical model due 
originally to Kamata & Nishimura (1958), and 
Greisen (1965), often referred to as the NKG ap- 
proximation. In the NKG model the longitudi- 
nal (along-track) development of the shower at a 
depth d is parameterized by its age s: 



s{d) 



3d/ Xo 



d/Xo + 2ln{E/Ecrit) 



(9) 



where Ecrit = 86 MeV for electrons in air, and 
Xq = 36.7 gm cm~^ is the electron radiation 
length in air. The evolution of the total num- 
ber of charged particles (virtually all electrons 
and positrons) with depth is then approximated 

by 



Ne = 



0.31exp[((i/Xo)(l - 1.51ns)] 
VHE/Ecrit) 



(10) 



Although this approach does not account for any 
of the large fluctuations that are possible in high 
energy air showers, it describes the average be- 
havior reasonably well. 

3.2. Transverse ionization density 

The transverse charged particle densities are 
described in a similar fashion, also parameterized 
by the age of the shower and the Moliere radius 
rm (Bourdeau 1980): 
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Fig. 1. — Electron ionization line density for 5 
showers of energies in the 10^^ to lO^i eV range. 
Such line densities are quite similar to those of 
radio meteors. 

and r is the gamma function, and s„j = 0.78 — 
0.21s. The Moliere radius for air is given by 

Vm = 2.12 X 10^ = 70 (^) m (13) 

Ecrit P \ P J 

where p is the density of the air at the altitude 
under consideration, and po is the sea level den- 
sity. The calculated density is in units of 
charged particles per unit area, passing through 
a plane transverse to the shower axis at the given 
depth. Since the shower thickness is typically of 
order several meters or less over its length, .^e 
can be approximately equated with the number 
density in a slab transverse to the shower. 

To convert these particle number densities to 
the resulting ionization density, we use the av- 
erage value for electron energy loss given by 
Ecrit I = 2.343 MeV gm"^ cm-^ divided by 
the mean energy per ion-pair for air Eion = 33.8 
eV (Segre 1977) which accounts for the ineffi- 
ciency of ion-pair production. This procedure is 
similar to that used by the Fly's Eye (Baltru- 



saitas et al. 1985) in estimating shower parame- 
ters based on fluorescence yield. 

We note here that the NKG approach and sub- 
sequent modifications to it do not accurately ac- 
count for the electron distribution at large radii 
from the core, and will in general underestimate 
the electron density beyond the Moliere radius. 
However, since the ionization density at large 
radii is low compared to the core, there is very 
little coherent contribution from these electrons, 
in spite of their large integrated number. Our re- 
sults remain basically unchanged by this effect. 

3.3. Numerical results 

Fig. 1 shows curves of calculated ionization 
line densities for several EAS of different primary 
energies over the range of 10^^ — 10^^ eV. The 
showers are assumed to be propagating horizon- 
tally at an altitude of 10 km, and the shower pa- 
rameters are corrected for the lower air density at 
this altitude. The bottom axis shows the along- 
track distance corresponding to the depth shown 
along the top axis. As noted above, the line den- 
sities at these energies correspond to typical line 
densities of radio meteors, which are detected at 
heights of 80-120 km. However, the lateral dis- 
tributions will be quite different as noted above. 
We note also that the showers are physically very 
long, with the region of shower maximum extend- 
ing over many km of physical distance for most 
showers at these energies. 

In Fig. 2 we show the effective line density as 
a function of the radius within which the electron 
density is integrated, for the inner 40 m radius 
for the same set of showers shown in Fig. 1. Here 
it is evident that the effective line densities of the 
inner ~ 10 m radius core of the showers arc about 
a factor of 10 lower than those which include all 
the electrons produced in the shower. However, 
this area corresponds to only 0.25% of the area of 
the Moliere disk ~ 200 m at 10 km), It is thus 
evident that the cores of the showers are highly 
concentrated in ionization density with respect 
to the bulk of the shower disk area. 
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Fig. 2. — The effective line density of an EAS 
near its shower maximum as a function of the 
radius within which the electron density is inte- 
grated. It is evident that the effective line densi- 
ties within radii of 10 m or so arc about 10% of 
the total line densities shown in Fig. 1, and are 
still substantial. 

In Fig. 3, we show the lateral ionization den- 
sity distributions near shower maximum for the 
same four EAS presented in Fig. 1 (solid lines). 
Included also are curves of the effective plasma 
frequency corresponding to the density at each 
radial distance (dashed lines). Given the uncer- 
tainty in the accuracy of the analytical model 
for EAS development at very small core radii (cf. 
r < 20 cm) , the implied highest radar frequencies 
that will undergo total reflection are in the range 
of 10-50 MHz, with a strong dependence on the 
primary energy. 

It is evident from Figures 2 and 3 that the 
cores of the showers are likely to be most impor- 
tant in determining the radar response, since the 
electron density is much higher near the core. In 
Fig. 4, we plot the core electron density again 
near shower maximum for a lO^'' eV horizon- 
tal shower, at a variety of different altitudes to 
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Fig. 3. — Radial dependence of the ionization 
density for the same showers presented in the pre- 
vious figure. Also shown is the radial dependence 
of the plasma frequency for each case. 

display the effects of the different atmospheric 
density. To first order the shape of the core is 
preserved with altitude, but the overall electron 
density decreases at higher altitudes because the 
production rate of ionization depends on the den- 
sity of the air. 

3.4. Duration of the radar echo 

Before estimating the radar echo power, we 
first address the important topic of the expected 
lifetime of the free electrons in the ionization col- 
umn. If EAS detection by radar is to be proven 
practical, it is first necessary to first establish 
that the ionization column can be treated as 
quasi-static over the timescalc it might take to 
interrogate the shower with a radar pulse, either 
as a standalone system, or if a separate air shower 
detector provided the trigger. 

As an example of the issue with respect to 
triggered operation, consider a large air shower 
array with radius 30 km. For typical electron- 
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Fig. 4. — Radial profile of showers shown for the 
inner 40 m radii, plotted on a linear distance 
scale, for different shower altitudes, showing the 
effects of altitude for a set of horizontal showers. 

ics used in such systems the trigger forms within 
several tens of //s. Round trip transit time to an 
EAS within ~ 30 km requires another 200 /its. 
Thus in this case we require EAS electrons to re- 
main mostly free for of order 250 /is or longer, if 
an external trigger is to be used to allow inter- 
rogation of a detected shower. In general, if the 
free electron lifetime is Tg, and the trigger for- 
mation and pulse initiation time is Tp, then the 
maximum range for a triggered system is 

Rtr < C . (14) 

If a radar detection system is operated in a 
standalone mode, the free electron lifetime will 
determine how often the detection volume must 
be pulsed in order to have a reasonable efficiency 
at detecting the EAS. In the limit where the free 
electron lifetime becomes so short that it is com- 
parable to the pulse duration, then pulsed radar 
detection will become impractical. Other meth- 
ods which use continuous- wave (CW) techniques 



may still apply in this case; these will be dis- 
cussed briefly in a later section. 

There are a number of effects that cause the 
electron density to evolve with time: 

1. Diffusion of the electrons through the am- 
bient air, lowering their density and thus 
their net radar cross section; 

2. Electronic recombination with ions; 

3. Attachment of electrons to neutral molecules, 
or dissociative attachment to atoms; 

4. Collisional detachment of previously at- 
tached electrons, which frees electrons that 
have been attached to atoms; 

3.4- 1- Diffusion effects 

Although the conditions of the atmosphere at 
the 80-120 km heights that meteor echoes are 
observed at are quite different from those at EAS 
heights, it is useful to preface our analysis with a 
discussion of the effects that determine the time 
scale for meteor echoes. 

The duration of meteor radar echoes has been 
modelled and studied experimentally in detail 
for many decades (cf. Hanbury Brown & Lovell 
1957; Kaiser 1968; Kaiser et al. 1969; Jones 
& Jones 1990; Jones 1991). The power of a 
meteor radar echo is found in the underdense 
regime to decay exponentially with a time con- 
stant Trn = A^/(327r^I?j) where Di is the ambipo- 
lar (or ion neutral) diffusion coefficient. At the 
altitudes that meteors are observed with radar, 
Dj~l — lOm^s"^, and the typical decay times 
for underdense trails are thus several tens of ms 
for frequencies in the VHF regime. 

To estimate diffusion effects at the lower alti- 
tudes of EAS ionization columns, we note that 
Di oc Tu~^ where T is the kinetic temperature 
and Ui the collision frequency (cf. Buonsanto et 
al. 1997). Thus at 10 km altitude, the diffusion 
coefficient is much smaller, Dj ~ 5 cm^ s~^, due 
mainly to the much higher collision frequency 
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at lower altitudes. The implied time constant, 
assuming diffusion progresses in the same man- 
ner for EAS as for meteor ionization columns, is 
Teas — 60 s for A = 3 m. 

S-band radar echoes from lightning at alti- 
tudes of ~ 10 km measured by Williams et al. 
(1984) were frequently found to persist for sev- 
eral hundred milliseconds. For a wavelength of 
11 cm, the diffusion time constant at this alti- 
tude is 240 ms, which suggests that diffusion may 
be the dominant effect in the evolution of light- 
ning radar cross section, at least at microwave 
frequencies. However, lightning ionization densi- 
ties are estimated to be initially 10^^ — 10^^ cm~^, 
with plasma frequencies 10^ or more times higher 
than for EAS columns. Thus other effects rele- 
vant to the conditions of this much denser plasma 
may be important. 

3.4-2. Attachment & Recombination 

Since diffusion is less effective at dissipating 
the ionization column of an EAS than it is for 
a meteor, at least at VHF frequencies, we ex- 
pect that attachment and recombination may 
play a more significant role. The recombination 
rate for atmospheric electrons is described by (cf. 
Thomas 1971): 

^ = - aeul - PeUe (15) 

where N^, rXg are the initial and evolving elec- 
tron density, respectively, is the electron-ion 
recombination coefficient (in units of cm^ s^^), /3e 
and is the electron attachment rate (s~^). Here 
we have ignored terms in Thomas (1971) which 
relate to the presence of negative ions in the iono- 
sphere since we are considering EAS within the 
troposphere where such ions are not prevalent. 

The electronic recombination coefficient for 
hydrogenic atoms can be approximated by (Seaton 
1959) 

tte = 5.20 X 10"^^ 
x( 0.429 0.5 In Q-F 0.469 Q'^/^ ) cm^ s"^ 



where Q = 1.58 x 10^ /Tg and Tg is the electron ki- 
netic temperature, with Tg ~ 10'^ — 10^ K typical 
for ionization in this case. For air in the tropo- 
sphere we thus expect cte ~ 10~^^ to 10~^^ cm^ 
s~^, and the implied recombination time for the 
highest electron densities we have considered is 
several minutes or more. Thus it is highly prob- 
able that the electrons will attach before recom- 
bining. ^ 

Molecular oxygen is electronegative, forming 
ions with free electrons primarily through disso- 
ciative attachment, via the process 

e" + O2 ^ O" + O (16) 

where the resulting oxygen ion and atom carry 
off the excess kinetic energy of the electron, and 
the electron affinity of O2 is of order 0.5 eV. 
Molecular nitrogen, in contrast, does not form 
attachments with electrons but can be more eas- 
ily collisionally excited, and may play a role in 
mediating detachment processes in air. 

Several studies of electron attachment in air 
using electron swarm techniques (see Gallagher 
et al. 1983 for a compilation) have indicated a 
much lower attachment rate than that which is 
expected from using the attachment coefficient 
for pure molecular oxygen reduced by the molar 
ratio in air (Moruzzi & Price 1974 &: references 
therein). Measurements of the attachment co- 
efficient are very sensitive to contamination by 
carbon dioxide, and CO2 was not removed from 
the samples in many early measurements. Al- 
though CO2 can be present in significant concen- 
trations in the near sea level atmosphere, its aver- 
age concentration in the atmosphere as a whole is 
0.033%. Thus electron attachment rates for dry, 
high-altitude air must be determined with CO2 
largely removed from typical air samples. 

^Wc note that values for the total recombination coeffi- 
cient of Q = 2 X 10~^ cm^ are found in the D region 
of the ionosphere (Thomas 1971). Such a high value is 
likely to be due to the large variety of ion species present, 
which increases the total recombination cross section per 
electron. 
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Moruzzi & Price (1974), after carefully remov- 
ing all CO2 from their samples, confirmed ear- 
lier measurements in being unable to detect any 
attachment in pure dry air. This result was at- 
tributed to a competing process: rapid collisional 
detachment from N2 interactions, yielding an ef- 
fective attachment rate that was less than 10% of 
the ionization rate and thus too small too mea- 
sure. Their limit on the effective attachment co- 
efficient [r]* /Njn) for pure air is 

{rf/Nm) < 3 X 10-^° cm^. (17) 

where Ny^ is the number density of the attach- 
ing molecular species. The attachment rate (3^ 
is related to the effective attachment coefficient 
through 

(3e = {t}* /N^{h))UeNra{h) (18) 

where Ue is the electron drift velocity. For a typ- 
ical electron drift velocity of tig = 2 x 10^ cm~^ 
s"-*^, the limit on r] gives f3e < 5.4x 10^ s~^, imply- 
ing that attachment dominates completely over 
recombination. 

Figure 5 shows the behavior of the electron 
density for three different altitudes for a 10^^ eV 
shower core region. Here we have plotted the 
electron density time evolution for both the case 
where the attachment rate is equal to the limit 
of Moruzzi & Price (1974), and one-tenth of the 
limit (dashed lines). It is evident that the ex- 
pected electron lifetimes are at least 10-20 /xs, 
increasing slowly at the higher altitudes. 

Although this estimate of the electron attach- 
ment lifetime is important from the point-of-view 
of evaluating whether there is sufficient time for 
triggered- interrogation of an EAS with radar, it 
is also important to establish an upper limit on 
the electron lifetime, both for purposes of evalu- 
ating the possibility of multiple-pulse interroga- 
tion of a shower, and for considerations of radar 
clutter from echoes of multiple showers for a 
radar system with very large range. 

The lightning radar measurements noted above 
(Williams et al. 1984 & references therein) of- 
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Fig. 5. — The evolution of the electron number 
density is shown for three 10^^ eV EAS, at differ- 
ent altitudes. The solid lines correspond to the 
evolution for the present upper limit on the at- 
tachment rate, and the dashed lines to one-tenth 
of the upper limit. 

ten found radar echo persistence of several hun- 
dred ms using repetitive radar pulses. However, 
the plasmas produced by lightning are extremely 
overdense and hot relative to those we are consid- 
ering, and thermal detachment may play a signif- 
icant role. In addition, the large amount of water 
vapor present in typical lightning radar measure- 
ments complicates any comparison with EAS in 
dry air. Thus comparison of radar echoes from 
lightning with those from EAS ionization must 
be done with caution. 

An rough upper limit for the free electron life- 
time can be obtained by considering the equilib- 
rium level of free electron content of the meso- 
sphere in the region of maximum cosmic ray 
ionization, peaking at about 12-14 km altitude. 
Here the dominant ionization is caused by low 
energy (several hundred GeV and above) cosmic 
ray air showers, which produce of order 30 ion 
pairs per cm^ per second (Rossi 1966; Gregory & 
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Clay 1982), roughly 15 times the sea level rate of 
ion production. 

The DC conductivity of the atmosphere at 12 
km is c7o ~ 4 X 10"^^ Q-^m'^ (Dolezalek 1982). 
The cosmic ray ionization at this altitude induces 
a tenuous plasma, with DC conductivity 



(To = 47reo — - 



(19) 



where eo is the dielectric permittivity constant, 
and the electron-neutral molecule collision fre- 
quency Ven can be written as 



4 X 10' 



-10 



% Nrr 



(20) 



where is the molecule number density (cm~^) 
at the altitude considered. The electron tem- 
perature Te is uncertain but likely to be in the 
range of 300-2000 K. At 12 km we thus expect 
z/ei ^ 4 X 10^° s^^ for Te = 300 K. 

Inverting equation 19 for the plasma frequency, 
which in turn depends on the equilibrium elec- 
tron number density (equation 1), we find 



Tie C± 35.5 Ven Cm 



(21) 



For the values determined above we have ng < 
0.57 cm~^ which implies a mean lifetime of less 
than 20 ms. The actual electron lifetime will be 
less than this since the detailed DC conductiv- 
ity also includes contributions from ionic terms 
which we have neglected here. There will also be 
an altitude dependence on this upper limit since 
it depends on the collision frequency. In sum- 
mary, we find that the available time for EAS 
interrogation by radar is bounded at 10 km alti- 
tude approximately by 



20 IIS < Te < 20 ms 



(22) 



Thus, for the shortest values of Te, triggered in- 
terrogation of air showers appears to be imprac- 
tical except for ranges of a few km or less; for 
the largest possible values, ranges of hundreds to 
even a thousand km or more are not excluded 



for higher altitude showers. In the latter case, 
no existing system can yet trigger on showers at 
these great distances, but planned space-based 
systems such as OWL/ Air Watch will have such 
capability. 

The ~ 20 ms upper limit to the electron life- 
time shows also that radar clutter from echoes of 
multiple showers should be completely negligible. 
For example, even assuming an energy thresh- 
old as low as 10^^ eV could be attained out to 
10 km, the shower rate at this energy over this 
area is still only several events per second. Radar 
systems also generally allow range-gated triggers, 
and this can facilitate rejecting any nearby clut- 
ter from low-energy showers within a few km of 
the transmitter. 

4. Predicted radar return power from EAS 

Having established that the duration of the 
EAS radar targets is likely to be sufficient for 
their detection, and introduced the radar cross 
sections expected, we now turn to estimation of 
the return power of the echo. 

Radar return power Pr is described in terms of 
a model where the radiation is emitted from an 
antenna with peak transmitted power P( and di- 
rectivity gain G = 47:^2^^, where Q.a is the solid 
angle of the main beam of the antenna. The ra- 
diation is assumed to then scatter from objects 
in the antenna beam and be re-radiated isotrop- 
ically in the frame of the scatterer, producing a 
R~'^ dependence in the returned power as a func- 
tion of range R. Deviations from isotropic scat- 
tering are thus absorbed into the effective radar 
backscatter cross-section ai,, which can be larger 
or smaller than the physical cross-section of the 
object. In addition any real transmitting and re- 
ceiving system will have less than unity efficiency, 
which we designate here as rj. The radar equation 
under these conditions is (cf. Skolnik 1990) 



Pt 



Here we are assuming that the transmitting and 
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receiving antennas are identical, and we are ne- 
glecting for the moment any polarization effects 
or losses in the medium. ^ 

Given equation 23, the problem of determin- 
ing the detectability of EAS-initiated ionization 
columns reduces to that of determining the ef- 
fective RCS (Tfe for a given choice of operating 
radar frequency, and the noise power of the spe- 
cific radar in use. The noise power is given by 
Pjv = kTsysAf, where Tgys is the system noise 
temperature, k is Boltzmann's constant, and A/ 
the effective receiving bandwidth, assumed here 
to be matched to the transmitting bandwidth. 

Pulse compression radar. Almost all mod- 
ern pulsed radar systems now use what is known 
as pulse compression, a method which allows the 
receiver bandwidth, and thus the noise power, to 
be minimized (cf. Wehner 1987). Pulse compres- 
sion is typically implemented by effectively dis- 
persing a band-limited pulse with an initial band- 
width A/o = (Afo)~^ through a filter, trans- 
forming it into a frequency chirp which spans 
the original bandwidth, but now has a duration 
At 3> AtQ. The receiver then uses an inverse 
filter to de-disperse the received pulse. The ef- 
fective bandwidth is thus A/ = {At)~^, but the 
range resolution AR cx Ato of the full bandwidth 
is recovered. 

CW radar. Continuous wave radar systems 
which transmit at 100% duty cycle, can derive 
range information by using frequency- or phase- 
encoded waveforms to uniquely tag the detected 
echoes to a particular time segment of the trans- 
mitted waveform. Processing of the received sig- 

*This equation is also derived strictly under conditions 
where the received radiation is in the far-field, that is, 
where R > 2D^ /X where D is the largest projected dimen- 
sion of the scattering target. In our case, this is not gen- 
erally satisfied, since the length of the ionization columns 
can be tens of km. However, we have already accounted 
for these Fresnel zone effects by limiting our analysis to 
the first Fresnel zone as noted above. 



nal typically involves mixing with a reference sig- 
nal, which transforms different ranges to different 
frequencies, followed by detection in a filterbank. 
The effective bandwidth is then determined by 
the bandwidth of the frequency channels used. 

For simplicity we will assume a pulsed radar 
system here and in what follows. 

Combining the noise power equation above 
with equation 23, the signal-to-noise ratio (SNR) 
of the received power is 

^ = "''^^^ kT^f ■ ^''^ 

Evaluating equation 24 for a nominal choice 
of parameters gives the SNR per received radar 
pulse per square meter of RCS: 

^ = 3 3^^^ (^] (^y 

N ' VI \1 kWj VO.iy \10j 

\3mJ VlO'^m; VIO^K; \WhsJ' 

(25) 

We note that the reference values chosen here 
represent a modest radar system; in particular 
the peak power of 1 kW is easily attained by cur- 
rent standards, and the directivity G ~ 10 (beam 
size ~ 1 sr) represents a relatively low-gain an- 
tenna. The system temperature Tgys = 1000 K 
is realistic for 100 MHz (A = 3 m) however, since 
the brightness temperature of the sky is quite 
high at these frequencies. 

Radar systems also routinely use repetitive 
pulsing to increase SNR, which then grows roughly 

1/2 

as Np where Np is the number of pulses that 
are averaged. In the case of an EAS, the number 
of pulse repetitions is limited by the free electron 
lifetime, with large uncertainty. For the short- 
est possible electron lifetimes, it is unlikely that 
repetitive pulsing can provide much improvement 
in SNR, but for lifetimes of order 1 ms, repetitive 
pulsing could provide substantial improvements 
in SNR. 
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Fig. 6. — Radar cross section over a range of an- 
gles centered on normal incidence (at 90°) for an 
air shower with Eq = 10^° eV, at a distance of 
10 km, at radar wavelengths 30, 15, and 10 m. 
Results are based on the thin wire approxima- 
tion. The polarization of the incident wave has 
been taken to be 1 radian with respect to the 
shower axis; thus these RCS values should apply 
on average to random relative orientations. 

We now consider in detail the expected RCS 
for EAS ionization columns in the various density 
regimes that were introduced above. 

4.1. Low- frequency /overdense case 

At radar frequencies below the plasma fre- 
quency of the ionized core of the EAS, the RCS is 
comparable to that of a metal cylinder, as noted 
above, with a radius of Tc and a length equal to 
that of the Frcsncl zone. Based on the analyti- 
cal shower model used above and the results in 
Fig. 2, the critical radius at a given wavelength 
and EAS cascade energy Eq can be empirically 



Fig. 7. — Normalized power loss factor due to 
loss of coherence for two horizontal E = lO^'^ eV 
showers at heights of 5 and 10 km. 
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It is evident from this result that for frequencies 
above 10 MHz, and for all but the highest EAS 

energies considered here, Tc <C A and we can use 
the thin wire approximation discussed above. 

In Fig. 6 we show results for radar cross sec- 
tions of the first Fresnel zone for a shower of 
Eo = 10^0 eV, for three different frequencies, for 
a range of angles centered on normal incidence, 
and a polarization angle of 1 rad with respect 
to the shower axis. The cross section at normal 
incidence is greatly enhanced due to the specu- 
lar reflection at that angle. We can empirically 
express the maximum cross section for this case 
as 



af{W km) = 2.2 x 10^ 
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m. (27) 
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It is clear from this analysis^ that overdense 
showers will show very large cross sections when 
the reflection is nearly specular. However, as seen 
in Fig. 6, the cross sections decrease rapidly 
for non-specular reflections. This behavior will 
be seen to contrast sharply with the case of un- 
derdense scattering, to be treated below. Any 
consideration of the detectability of a particular 
shower must account for both contributions to 
the cross section. 

In the analysis above we have also ignored the 
additional contribution due to the partially co- 
herent scattering of the much larger volume of 
underdense plasma outside of Vc- This region 
could have refractive effects on the incoming ra- 
diation. In an appendix we treat this issue in 
some detail and show that it produces a negli- 
gible effect for the showers under consideration 
here. 

4.2. High- frequency /underdense case 

When either the frequency is high enough or 
the ionization density low enough that there is ef- 
fectively no region of the column where the crit- 
ical density obtains, the radar return is due to 
Thomson scattering of the free electrons in the 
column, modulated by the phase factor associ- 
ated with the physical extent of the ionized re- 
gion as described in equation 8 above. 

We have made numerical estimates of the ef- 
fective cross section for a restricted set of cases 
here, with the methods described in more detail 
in an appendix. Here we summarize the results. 

4-2.1. Normal incidence 

For the case of radar observations at normal 
incidence to the shower, the near constancy of 

®Here the linear increase of the cross-section with shower 
impact parameter R is due to the fact that the size of 
the first Fresnel zone grows with distance. In fact, al- 
though we have used the length of the Presnel zone Lf 
as the characteristic length of the "thin wire" portion of 
the shower, the actual length of this region is likely to be 
much longer than Lp. 



the ionization density along the track allows us to 
estimate the frequency dependence of the phase 
factor via a Fast Fourier transform (FFT) of 
the two dimensional ionization distribution in a 
plane transverse to the track. We have taken 
advantage of this feature to make numerical es- 
timates of the frequency-dependent behavior of 
the phase factors at normal incidence. 

Note that the requirement of normal incidence 
for the radar rays interrogating the shower is not 

as restrictive as it might seem. For example, con- 
sider a horizontal shower whose maximum is at a 
range of Ri from the radar transmitter /receiver 
system, and whose axis is inclined an angle 6 with 
respect to the ray that joins the radar source and 
the shower maximum. Under these conditions, 
the radar system can still interrogate the shower 
along a ray that intersects the shower at normal 
incidence at a distance Ri cos 6 from the shower 
maximum. 

Thus a horizontal shower at 10 km altitude, 
with 15 km range to its maximum, and inclined 
at ^ = 60° will have a ray that intersects it at 
normal incidence at a longitudinal distance along 
the shower of 7.5 km from shower maximum, at a 
range of 13 km, where the ionization line density 
is still ~ 60% of its maximum value. Thus the 
corresponding acceptance solid angle over which 
these results apply is not negligible, if the beam 
of the radar system is not too restrictive in solid 
angle. 

Fig. 7 shows our estimates of the frequency 
dependence of the phase factor for two horizon- 
tal showers at altitudes of 5 and 10 km, for 
£"0 = lO^'^ eV. The phase factor depends strongly 
on both frequency and shower altitude. The fre- 
quency dependence is due to loss of coherence as 
noted above, and the altitude dependence arises 
from the change in Moliere radius to smaller val- 
ues at lower altitudes. We have also investigated 
the behavior at other energies in this regime and 
found these results to be insensitive to energy 
over the range 10^^ — 10^^ eV. For the two cases 
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shown we have included fitted power laws 



h=10 km, Eq = 10^° eV 



$(/;5km) = 0.114/foi-i5 

$(/; 10 km) = 1.13 X lO'^ /-^-^^ . 

Here /lo is the radar frequency in units of 10 
MHz. These relations are individually accurate 
to about 20% over the range from 10 MHz to 0.3 
GHz. 

Given the power loss associated with the phase 
factor, the total cross section is then estimated by 
summing the cross sections of the individual vol- 
ume elements over one Fresnel zone of the track, 
by analogy to the overdense case, and multiply- 
ing by the overall phase factor for the radar fre- 
quency and shower altitude. 

Evaluating this modified cross section for its 
energy and frequency dependence, we derive an 
empirical expression for the case of a horizontal 
EAS at altitude of 10 km 

.-(10 km) = 175 (^) 

and similarly for an altitude of 5 km: 

..-(5 to) = 1400(3^)"' " 

'<(i()^iv)''°(iol^) 

These relations should be valid out to ranges 
of several tens of km. At the highest energies and 
lowest frequencies the underdense case may not 
obtain, but these estimates will then provide a 
lower limit on the RCS. 

4-2.2. Oblique angles 

For shower angles away from normal inci- 
dence, it is necessary to consider all of the in- 
dividual scattering contributions of the electrons 
over the entire volume of the shower, and to treat 
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Fig. 8. — Radar cross section over a range of an- 
gles centered on normal incidence (at 90°) for a 
horizontal air shower with Eq = 10^*^ eV, at an 
altitude of 10 km, at radar wavelengths 30, 10, 
and 6 m. Results are based on numerically inte- 
grating the contributions of all of the individual 
volume elements of the shower. 

the pulse echo behavior explicitly. We discuss 
this in detail in an appendix. We have numer- 
ically estimated the integrated RCS for several 
wavelengths at a distance of 10 km to display 
some of the characteristics of the angular depen- 
dence. 

The results of this are shown in Figures 8 and 
9. We have plotted the total RCS as a func- 
tion of angle over about 1 radian around normal 
incidence, for a horizontal shower at heights 10 
and 5 km and energy 10^*^ eV, and several dif- 
ferent radar wavelengths. The dependence on 
angle is quite different from the overdense case 
where specular reflection is important; here, the 
cross section has an approximately Gaussian pro- 
file with angle, with a FWHM of order 30°. The 
fairly steep wavelength dependence is also inde- 
pendent of angle over this range. Thus the empir- 
ical relations derived in the previous section can 



16 



h=5 km, = 10^° eV 




angle (degrees) 



Fig. 9. — Similar to previous figure, with shower 
altitude of 5 km. 

be used with appropriate factors for the angle of 
incidence shown here. 

We caution here that the angular dependence 
of the RCS at distances that are much greater 
than the length of the shower is probably not 
accurately reflected by these results. We have 
not yet derived a general analytical treatment of 
the behavior for large distances, but we note that 
the methods described in the appendix allow for 
numerical estimates to be made at any distance. 

4.3. Caveats 

We note that, due to self-imposed limitations 
of scope, a number of effects have not been 
treated in our analysis. We list several of these 
here. 

Geomagnetic effects. We have neglected ge- 
omagnetic charge separation in the shower de- 
velopment. Although this is an important effect 

in determining the detailed shower structure, it 
depends greatly on the shower direction with re- 
spect to the magnetic field. Its effect on the RCS 



would also depend on the angle of observation 
and the treatment of this complex interaction is 
beyond our scope at present. 

Landau-Pomeranchuk-Migdal (LPM) effect 

We have neglected the LPM effect in the showers 
considered here. This effect, which retards the 
initial growth of the shower, is important at the 
energies discussed and should be accounted for 
in a more detailed treatment. However, it affects 
primarily shower which interact initially through 
the electromagnetic rather than the hadronic 
channel, and thus there are still many EAS at 
these energies for which the approximation used 
is still accurate. 

Secondary Fresnel zone echoes. For sim- 
plicity we have only treated the echo from the 
primary Fresnel zone here. In fact, for the high 
SNR case, many distinct echoes may be detected 
from successive Fresnel zones, and these will cer- 
tainly add significantly to the understanding of 
any detected shower. On the other hand, for 
undcrdense showers that are not observed near 
normal incidence, the time spread of the return 
will also tend to decrease the instantaneous SNR. 
Such effects are important in understanding the 
response of a radar system to the wide variety of 
possible EAS directions and relative angles and 
will need to be explored further in any real ap- 
plication. 

Near vertical EAS. For convenience we have 
analyzed a horizontal air shower, but such show- 
ers are of general interest, since such highly- 
inclined showers imply very deep initial interac- 
tions and are thus a possible signature of neu- 
trino primaries. For more typical showers with 
small zenith angles, the Moliere radius changes 
significantly during the shower development and 
complicated the analysis considerably. General- 
ization of our analysis to EAS at all angles is best 
done with a full numerical simulation. Prelimi- 
nary results of such a simulation (D. Bergman 
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2000, pers. comm.) do not show significant dif- 
ferences in the ionization profiles, and still show 
a tight core in the radial ionization density. 

Pulse smearing by the shower. Note that 
our results for the under dense case above apply in 
a general way only to non-pulsed radar, since the 
superposition implied by the Fourier transforms 
assumes that each Fourier component extends 
spatially out to dimensions that are much larger 
than the shower column. For pulses whose com- 
pressed time scale is physically shorter than the 
time-projected thickness of the ionization column 
of a shower, wc must explicitly account for the 
behavior of the individual pulse echoes from dif- 
ferent portions of the shower. However, the effec- 
tive RCS corresponding to the integrated power 
of the echo is still described by our analysis to 
first order. 

A general consideration of the complications 
arising from the pulsed nature of the radar are 
beyond our scope at present but must be included 
in any detailed estimate of the pulsed radar echo 
behavior. In the appendix we describe a method 
developed to deal with pulsed-radar reflections, 
and we utilize this method in some later results. 

Plasma resonance effects. We have treated 
the undcrdcnsc case in a simplified way which 
neglects the possibility of plasma resonance ef- 
fects. Such effects are in fact seen in meteor 
echoes (Poultcr & Baggaley 1977) and can play 
an important role in enhancing the RCS, and in 
producing more complex time structure in the 
radar return. 

5. Discussion 

In the previous sections we have presented 
evidence that FAS ionization columns are well 
within the range of VHF radar detection for cas- 
cade energies above ~ 10^^ eV. Now we turn to a 
discussion of the applicability of this approach to 
present efforts at detection and characterization 



of such FAS. To establish a basis for comparison, 
we first outline the capabilities of fluorescence de- 
tector systems, to which FAS radar detection, if 
possible, would be most similar. 

Fluorescence detector capabilities. The most 
capable present fluorescence detector system is 
the High Resolution Fly's Fye system in Utah 
(Abu-Zayyad et al. 1999) operating for several 
years now in "monocular" mode, and now begin- 
ning to come on line with stereo-mode operation, 
using a pair of telescope clusters separated by 
12.6 km. HiRes uses about 1 photomultiplier per 
square degree of focal-plane sky coverage, with 
2-3 m aperture low-resolution optical telescopes 
used as the light collectors. 

HiRes can detect events out to 20-25 km, al- 
though at this distance the effects of aerosols 
increase the uncertainty of the measurements. 
To achieve reasonable precision in the measure- 
ment of 6 shower parameters (three geometric 
and three associated with the energy and shower 
development) requires typically a minimum of 
about 20 photomultiplier pixel hits in monocular 
mode, giving at this level an energy resolution of 
~ 25% and precision of about 30 gm cm~^ in the 
position of shower maximum (R Sokolsky 2000, 
pers. comm.). Significant improvement in the 
precision of these measurements is expected once 
the stereo mode comes into full operation (Abu- 
Zayid et al. 1999). Observation efficiency for any 
Fluorescence detector is limited to about 10% 
since the observations can only be made on clear, 
moonless nights. The total acceptance aperture 
for HiRes is expected to be about 6000 km^ sr in 
the region of 10^*^ eV shower energy. 

5.1. Standalone EAS radar system. 

Here we consider what type of radar system 
would be required to produce EAS measurements 
that can be compared to those expected from 
HiRes. 

Because of the issues of interference rejection. 
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a standalone EAS radar system would almost 
certainly have to be constructed with multiple 
stations configured to both transmit and receive, 
with pulse encoding such that each station would 
be sensitive to pulses produced by the other sta- 
tions. Bistatic radar cross sections of EAS will 
vary only modestly from the monostatic RCS val- 
ues treated here. Global synchronization of the 
stations to several tens of ns or less is routine 
with the Global Positioning System, particularly 
with the disabling of Selective Availability (SA) 
of the high precision GPS codes. 

To optimize the signal-to-noise ratio of the de- 
tected echoes, we wish to operate at the low- 
est frequency that is practical from considera- 
tions of background noise, including ionospheric 
noise from daytime operation. Although lower 
frequency operation may be possible, we choose 
here to consider only frequencies above about 30 
MHz. In the HF-VHF regime the effective system 
noise temperature for a remote site (distant from 
urban noise sources) can be written as (Skolnik 
1986): 

which yields T^ys = 3600K at 30 MHz. 

Assuming that each station receives an echo 
from its own as well as all others in an ar- 
ray of M transmit/receive stations, there are 
M + M{M — l)/2 complex measurements for 
each combined pulse of the M stations. For 
six EAS parameters to be fitted, the minimum 
number of stations required is M = 3 (giving 
6 total geometric ranges and echo amplitudes), 
assuming that the measurements are indepen- 
dent. Given that there will be some covariance 
of the measurements, and requiring that the fit 
be well-constrained, at least 4 or 5 stations will 
be needed, giving 12-20 ranges and amplitudes. 
Additional closure-delay (or triangle) measure- 
ments number M(M — 1){M — 2)/6, and will 
provide an extra 4-10 independent constraints for 
the M = 4, M = 5 cases, respectively. However, 



these latter constraints apply only to the geom- 
etry, which is typically well-constrained in any 
case, and not the echo amplitude. Addition of 
dual-polarization capability, which is not uncom- 
mon in radar systems, would double the number 
of measurements. 

5.1.1. Radar data types. 

The measured data can be characterized as 
follows: 

Range data &: shov^er angle. The time delay 
to the leading edge of the echo at each station 
will provide the range from the transmitter to the 
receiver. The ranges can be measured typically 
to a precision of roughly the compressed pulse 
width divided by the SNR. For a SNR=10, and 
a compressed pulse width of 250 ns, the range 
resolution is below 10 m. 

A complete determination of the angular res- 
olution possible using radar is beyond our scope, 
but simple arguments can show that a set of 4 or 
5 stations should provide reasonable estimates of 
the track direction. A schematic view of the ge- 
ometry associated with the range measurements 
for a pair of stations is shown in Fig. 10. Here 
the stations 1 and 2 measure leading-edge ranges 
Ri and i?2 to separate portions of the shower 
at along-track positions si,S2 with respect to 
shower maximum. The stations are separated 
by a baseline B12 which is not generally parallel 
to the shower track. Let a vector n bisect the 
baseline and the shower segment from si to S2- 
The shower direction can then be characterized 
by angles a, which measures the rotation of the 
shower track around the axis n with respect to 
the projected baseline, and /? which measures the 
angle between the shower and n. 

It is evident from Fig. 10 that changes in j3 
will under most conditions produce first order 
changes in the range values, eg. 

|S2-Sl| 
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shower track 




station 1 



Fig. 10. — A schematic of the range constraints 

on shower angles provided by a pair of radar sta- 
tions 1 and 2, separated by basehne B12. See 
text for details. 

where aR is the mean range error and cr^ the 
resulting angular resolution in (3. Since for tracks 
that are primarily transverse to the baseline B12 , 
we have |s2 — si| ~ B12, then baselines of order 5 
km or more with range errors of ~ 10 m will give 
angular resolution in /? of order several mrad, or 
about 0.1°. 

Changes in a for a pair of stations under typi- 
cal expected detection geometry do not give first 
order effects in range, and are thus more prob- 
lematic. Because a change in a corresponds to 
a kind of "torsion" mode of the range contraints 
(using a rigid-rod analogy), the error in a for a 
two station solution is likely to generally behave 
as 

(Ja ~ 1 f (^2j 

\S2 - Sl\ 

which applies most directly to the case where 
a <C 1. For example, for i? = 20 km, As = 5 

km, and (T/j = 10 m, we have CTq ~ 7°, much 
lower precision than can be achieved in air flu- 



orescence detectors. We note however that the 
addition of more radar stations, with the inclu- 
sion of bistatic ranges and multiple polarization 
measurements should rapidly improve the preci- 
sion as more constraints are added, but further 
estimates are beyond our scope at present. 

Echo amplitude data. The echo amplitude is 
more difflcult to interpret. For under dense show- 
ers, the echo amplitude is directly proportional 
to the number of electrons present, and thus is 
similar to fluorescence detectors in terms of its 
ability to estimate the number of ions within a 
certain region of the track. However, since the 
radar echo may be spread out over a band of dif- 
ferent ranges, determination of the contributing 
portion of the shower for each portion of the re- 
turned echo at each receiving station is a clearly 
more complicated problem than that faced by flu- 
orescence detectors, which rely on imaging preci- 
sion to solve this problem. However, since the ge- 
ometric parameters of the shower track are easily 
separable via the range data, the resulting am- 
plitude of the leading edge of the return can be 
attributed unambiguously to the portion of the 
ionization column that is nearest to the receiving 
station. By thus reconstructing a profile of these 
prompt received amplitudes at their relative po- 
sitions along the geometric track, the profile of 
the shower development can be extracted. 

For overdense showers, the radar echos may 
show anomalously high strength compared to ex- 
pectations based on an underdense echo. In this 
case, interpretation of the data may rely on dis- 
tinguishing these two cases on an event-by-event 
basis. Estimates of the effective conductivity of 
the plasma can then still yield the electron con- 
tent of the shower. 

In either of these cases, the amplitude data, 
although not as straightforward to interpret at 
that of fiuorescence measurements, still provides 
the type of calorimetry of an EAS that is neces- 
sary to determine the shower energy. This is a 
feature that particle counter arrays cannot match 
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in detail, and it sets apart the fluorescence de- 
tection and potential radar detection techniques 
from other methods. 

5.1.2. Duty Cycle. 

Perhaps the most interesting feature of radar 
as a potential EAS detection approach is the pos- 
sibility of much higher duty cycle observations 
over detection volumes comparable to those of 
fluorescence detectors. In principle a radar sys- 
tem is not limited to night operation, and suffers 
very little attenuation from clouds or aerosols. 
Although operation in the presence of lightning 
within the detection volume might produce some 
spurious events, avoiding operation in periods of 
lightning would cost at most a few percent of the 
livetime. 

The largest uncertainty in efficiency is proba- 
bly due to the free electron lifetime, since if this is 
very short, it may be difficult to ensure adequate 
pulse repetition frequency to detect all showers. 
And although radar propagation is unaffected by 
aerosols and clouds, the free electron lifetime is 
not immune to variations in these components 
of the atmosphere, and in practice they may be 
important limitations to the efficiency. In addi- 
tion, the presence of the sun above the horizon 
will increase the system temperature somewhat, 
especially near solar maximum, and this will act 
to increase the detection threshold for daytime 
operation. 

5.2. SNR Budget for a standalone EAS 
radar system 

In Figure 11 we show a plan view of the geom- 
etry of a radar array using both monostatic and 
bistatic detection of an EAS, here shown with 4 
stations at 3-5 km spacing, providing 10 ranges 
and 10 pulse amplitudes. Here we describe the 
parameters required of each transmit/receive sta- 
tion to satisfy the needs of the array described. 

If we assume that we require an energy thresh- 
old of 10^^ eV out to a range of 20 km, the im- 




3-5 km 



Fig. 11. — An array geometry for a possible stan- 
dalone EAS radar detection system. Each of the 
4 stations has both transmit &; receive capabil- 
ity, and can distinguish the signals from the other 
stations as well. 

plied RCS for the underdense case is ~ 2 m^ 
at this energy. The Fresnel zone length at this 
range is about 250 m for a shower at 10 km al- 
titude, corresponding to about 10 g cm^^ along 
the track. Using a commercial VHF radar sys- 
tem, it is straightforward to achieve 60 kW peak 
power and pulse repetition rates of 10-50 kHz for 
10 IIS pulses. The sky brightness temperature 
will lead to a system temperature of Tgys = 3600 
K at 30 MHz as noted above. Since we wish to be 
sensitive to as large a volume as possible, we use 
an antenna with a broad beam (thus correspond- 
ing to a low gain). A suitable system is a verti- 
cal monopole or discone system (Balanis 1997) 
which provides good 27r azimuthal coverage from 
5° to about 50° elevation, giving an overall gain 
of ~ 3 with respect to an isotropic antenna. We 
also assume here a somewhat pessimistic overall 
transmit /receive efficiency of r/ = 0.05, and a sin- 
gle pulse echo (no pulse averaging) to account for 
possible rapid decay of the free electron popula- 
tion. 
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Table 1: SNR budget for one station of a possible 
EAS standalone radar system. 



PflrflTTlPt'PT' 


V3illlG 


±dBm 


Received Power 












Peak transmit power 


60 kW 


77.8 


Pulse duration 


10/is 




chirp bandwidth 


3 MHz 




Number of repetitions 


1 


0.0 


Antenna gain 


3 


9.54 


wavelength 


10 m 


20.0 


atat E = 10^^ eV 


3.8 m^ 


5.8 


range to EAS 


20 km 


-172.0 


xmit/rcv efficiency 


0.05 


-13.0 




5.04 X 10~* 


-33.0 


Received power 




-104.9 dBm 


Noise power 






Boltmann's constant 


1.38 X 10"^° 


-198.6 




mWK-^Hz"^ 




System temperature 


3650 K 


35.6 


effective bandwidth 


100 kHz 


50.0 


Noise power 




-113.0 dBm 


SNR 


6.4 


8.0 dB 



Using these assumptions we show the SNR 
budget, in dBm (decibels referenced to 1 milli- 
watt) in Table 1. This budget corresponds to a 
single monostatic or bistatic pulse measurement 
in the array. The resulting SNR of 6.4 provides 
a range resolution of 34 m at a 30 MHz operat- 
ing frequency, using a chirp bandwidth of 10% (3 
MHz in this case). 

We can also estimate the rate of shower de- 
tection as a function of energy. To do this we 
estimate that the acceptance solid angle of the 
system is '-^ 1 sr, based on the RCS within 10° of 
the peak values. This is significantly lower than 
that of a fluorescence detector, and reflects the 
rapid decrease in the RCS as the angle of inci- 
dence moves away from the normal to the shower 



axis. We also consider only the RCS for showers 
at 10 km altitude, although showers at lower alti- 
tudes will have significantly higher cross sections 
in general. The integral number of events seen 
above a given energy Ethr is given by 

Q{Ethr) - T^Rlas riobs T Q I{E > Ethr) (33) 

where riobs is the observing efficiency, T the ob- 
servation duration, f2 ~ 1 is the effective solid 
angle, and 

^(^>^*/^^)-(l^)"km---V-^ 

(34) 

is the integral cosmic ray spectrum at these en- 
ergies. Here the maximum range for detection is 
determined by inverting equation 4 above, which 
yields: 



V 10 y V 1020 eV 



p 



1 kW 



1/3 



^sys 



VlOOO K 



_r]_ 
0.1 

1/3 



1/3 



At 



G\2/3 

w) 

1/3 



(35) 



,10 /xs. 

This relation is approximately valid to ener- 
gies of order 10^*^ eV, where the detection range 
is ~ 80 km. Beyond this range, the minimum 
beam elevation is such that it exceeds the initial 
assumption of 10 km shower altitude, and the de- 
tection efficiency function of shower altitude 
must be more explicity addressed However, this 
exercise highlights the large potential detection 
area for a radar system 

In Figure 12, we plot the expected event rate 
per year for 80% observing efficiency, which is 
achievable with a fully operational system. Here 
the assumed minimum SNR = 6.0. Also plot- 
ted with respect to the right axis is the effective 
aperture, in km^ sr, for the same assumptions. 

It is evident from this analysis that EAS de- 
tection by radar certainly merits further inves- 
tigation as a high-duty cycle alternative to ex- 
isting EAS detection systems. Although it is 
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Fig. 12. — The expected event rate detected at 
a SNR=6 for the standalone EAS radar system 
described in the text and table 1. The effective 
aperture corresponding to this rate is also plot- 
ted. 

not a competitive technique at lower energies, it 
appears that this approach may have significant 
benefits for detection of EAS at the highest en- 
ergies. Although it may not be able to compete 
with fluorescence detection in its ultimate accu- 
racy in the estimation of shower parameters, a 
system which utilizes a trigger provided by an 
existing particle detector array such as the Auger 
Observatory could potentially gain valuable com- 
plementary information on daytime events where 
fluorescence information cannot be obtained. 

6. Conclusions 

We have demonstrated that, using standard 
models for the average behavior of extensive 
air shower development, the resulting ionization 
appears straightforward to detect using radar 
techniques in the VHF frequency range (30-100 
MHz), for primary energies greater than about 
10^^ eV. We estimate that a relatively modest 



ground-based system could provide information 
on EAS that is comparable to that available from 
other techniques, and may do so with a much 
higher duty cycle and potentially a larger effec- 
tive aperture. The approach rcqiiircs no new 
technology, making use of standard radar equip- 
ment and detection techniques, and may thus 
also provide a low-cost alternative to existing 
methods. 
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a very valuable critique. We thank Douglas 
Bergman, George Resch, Pierre Sokolsky, Paul 
Sommers, David Saltzberg, and Trevor Weekes 
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A. Appendix 1: Numerical Evaluation of 
RCS for the underdense case 

To estimate the RCS for the underdense regime 
of an EAS ionization column, we need to evalu- 
ate the individual contributions of the scattering 
electrons and the phase of the resulting scattered 
fields. As noted in section 2 above, under cer- 
tain conditions this is equivalent to determining 
the Fourier transform of the electron distribu- 
tion. Here we describe the basis for this con- 
clusion and provide details of the methods used 
for numerical evaluation of the RCS under these 
circumstances. 

A.l. Geometry of the problem 

The geometry for this evaluation is shown in 
Fig. A13. Here we let R represent the vector 
distance from the origin O (at the radar trans- 
mitter/receiver) to the position of shower max- 
imum. We choose the x — y plane to be coin- 
cident with the plane containing the axis of the 
shower and R. Within the ionization column it- 
self, wc use cylindrical coordinates p, (t), z\ to cal- 
culate the proper ionization density at position 
P with respect to the position of shower maxi- 
mum; r' = (p, 0, z\) is the vector separation of P 
from shower maximum. The vector joining the 
origin to P is denoted as r. The angle between 
the shower axis and R is denoted d. 

Given these definitions, we find 

|r p = -I- 2R{p cos sin 6* - zx cos Q) 

(Al) 

which expresses the distance to an arbitrary 
point in the ionization column as a function of 
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shower track 

Fig. A13. — The geometry for numerical eval- 
uation of the summed radar cross section of an 
ionization column for the case of oblique angle 
of incidence. The radar transmitter /receiver is 
at the origin O, and the angle of the track is Q 
with respect to the direction from the origin to 
shower maximum R. The vector to the position 
P in the ionization column that is being evalu- 
ated is r, and the vector separation of P from 
the position of shower maximum is r'. is the 
azimuthal angle of P in cylindrical coordinates 
with respect to the shower axis. 

the distance to shower maximum, the local cylin- 
drical coordinates of the point, and the angle of 
obliquity Q of the shower. 

A.2. Scattered field 

The electric field that is incident at P can be 
written as 

Vine = UtrGtre'^'^ (A2) 

where Utr is the magnitude of the transmitted 
field, Gtr is the transmitter gain factor (as de- 
fined in the main text), k = 2tt/X, and r is de- 
fined in Fig. A13. Assuming that k and r are 
parallel (that is, neglecting the effective dielec- 
tric behavior of the underdense plasma in the col- 
umn), the received field from a small volume of 
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electrons UedV which scatter coherently is then 
given by 

dUrcv = UincV^nedVGsce^^— — 7- (A3) 

\^sc\ 

where Ggc is the effective gain of the electron 
scattering process [Gsc = 1 for isotropic scattering),^ 
and k^c, t^sc are the wavevector and radial vec- 
tor of the scattered wave. In this case we have 
kgc = — k and Vgc = — r, and thus 

dUrcv = UtrGtrGsc^f^riedVe^'^^ (A4) 

which displays the expected radial dependence 
of the field (equivalent to a r~'^ dependence in 
received power), and also the additional factor 
of two in the spatial phase which is due to the 
two-way trip of the radiation. 

If we assume for the moment that the time 
variation of the field is steady-state, we can drop 
the time dependence, and integrate over the vol- 
ume of the shower to yield the total received field 
strength, assuming that the transmitted beam is 
larger than the dimensions of the shower. If we 
also assume that the distance from the shower is 
large enough that |r| k, |R|, then we have 

R^Urcv = UtrGtrGscV^ [ ne{v)e''^-'^ d^Y 

Jv 

(A5) 

which shows (taking q = 2k) that the received 
field is proportional to the Fourier transform of 
the electron density distribution. 

A. 3. Normal incidence 

Under conditions where the pulse is incident 
on the ionization column at normal incidence, 
and if we assume that the ionization column is 
essentially constant along its longitudinal length. 



'For Thomson scattering, Gsc — 1.6 corresponding to elec- 
tric dipole emission. However, this response is comph- 
cated in any real plasma by a number of effects, and we 
have assumed Gsc = 1 for all cases considered here. 



the received field becomes proportional to the 
two-dimensional Fourier transform of the cross- 
sectional electron density, times the length of the 
Presnel zone L p (defined as the length over which 
the scattered fields are in phase): 

R^Urcv\e=o = UtrLpGtrGscV^ / ne{r)e''^-'' dh . 

J A 

(A6) 

We have numerically evaluated this using Fast 

Fourier transform techniques, since the area A 
over which this is done is of order A = Trr^ 
and the number of Xmin/^ cells in the grid is 
64^/A^j„, where we have assumed Xmin ~ -5 m, 
giving a can 
be normalized to provide a frequency-dependent 
phase factor, which we have plotted in Fig. 7 
above. 

A. 4. Pulsed radar &; the time-projected 
RCS 

In the previous analysis we have assumed a 
steady-state incident field (which corresponds to 
carrier-wave or CW radar). However, in our 
analysis of section 4, we assumed that the radar 
is pulsed, and in fact pulsed radar is preferable 
to achieve the kind of time resolution required 
in this application. For the ovcrdcnsc case, the 
analysis presented in section 4 is adequate, since 
we are considering refiections from a localized 
source at the critical radius of the ionization col- 
umn. Also, for normal incidence, the reflections 
from the first Frcsnel zone arrive closely in phase 
and can be treated without explicitly accounting 
for the pulse echo behavior. 

When we allow for oblique shower angles, how- 
ever, the concept of the Fresnel zone loses its use- 
fulness, and the received echo comes from the 
entire volume of the shower and is thus signifi- 
cantly spread out in time. We require a method 
of treating the behavior of a discrete pulse under 
these conditions. 

The approach we have followed is intuitively 
simple: we first produce a time-projected profile 
of the electron density; that is, we calculate the 
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ionization column 




radar transmitter/ 
receiver 

Fig. A14. — Illustration of the geometry that 
gives rise to the concept of the time-projected 
radar cross section. The ionization column can 
be divided into radial shells such that all of the 
electrons that scatter the incoming radar are in 
phase (to within a radian) within the shell. 

round-trip time associated with the global phase 
of each of the contributing differential volume el- 
ements of the scattering electrons 



Here ti is the round-trip time for the ith volume 
element at location Fj. If we consider all of the 
volume elements for which ti = tj within a spec- 
ified interval At < oj~^, then we define a shell at 
distance r = t/2c with a thickness of 1 radian of 
phase at the radar angular frequency lo. This is 
depicted in Fig. A14. 

Summing all of the contributing elements in 
this shell gives the time-projected square root of 
the total cross section at that phase: 

Y{t)At = ^ / ne{r)(fr . (A8) 

Jr-t/4c 



Once Y{t) is determined for a given radar 
wavelength, shower angle and energy, and dis- 
tance, then the pulse response of the ionization 
column is determined by convolving the electric 
field of the pulse profile with Y{t), and squaring 
the result. If the pulse is defined to have unit 
amplitude, then the resulting convolution can be 
expressed as an effective radar cross section for 
use in equation 4. 

For example, if the pulse is assumed to be a 
5- function in time, the resulting convolution T(i) 
is 

/oo 
Y{t')S{t - t')dt' = Y{t) (A9) 
-oo 

and the received power is then proportional to 
For other pulse profiles we numerically 
evaluate the convolution explicitly, and take the 
maximum of the resulting profile to be the effec- 
tive radar cross section. 

In Fig. A15 we show an example of the square 
of Y{t), giving the time-projected cross section 
calculated according to the manner described 
here. In each case the time bins correspond to 
a radian of angular phase according to the wave- 
length indicated, and the value for each bin is 
estimated according to equation A8. The pro- 
file in this case shows a characteristic rapid rise 
on its leading edge, reflecting the rise in cross- 
sectional ionization density of the nearest portion 
of the shower. The more slowly falling trailing 
edge is due to the contributions from portions of 
the shower further away. 

We note that Y{t) defined in this way is analo- 
gous to a Green's function for the time response 
of the system. It can be used in this manner 
to evaluate arbitrary radar pulse shapes, for a 
given EAS energy and direction. For extremely 
broad-band radar systems, the effective RCS val- 
ues will begin to approach those of Y{t). For 
more narrow-band radar pulses, the convolution 
of the field oscillations with Y{t) will act to par- 
tially cancel and the effective RCS will decrease. 
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Time-projected cross— section 
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Fig. A15. — The square of the function Y{t), 
which gives the time-projected total radar cross 
section, for an ionization column produced by a 
shower with Eq = 3 x 10^° eV, for A = 10 m, at 
a distance of 20 km. The time axis is windowed 
around the arrival time of the round trip travel 
time. Y{t) may be thought of as a kind of Green's 
function time response of the ionization column 
to a sharp radar pulse, with a magnitude that is 
proportional to the square root of the RCS at a 
given pulse phase. 

B. Appendix 2: Refractive effects of the 
plasma halo 

The refractive index of a plasma is given by 
Spitzer (1962) in the quasilongitudinal approxi- 
mation by 



n 



l-(^ 



1/2 



(Bl) 



where we have assumed that the frequency of the 
radiation is well above the electron gyrofrequency 
in the earth's magnetic field (~ 840 kHz). In 
this case the phase delay of the radiation after 
propagating a distance I through the plasma is 



^iven as (cf. Sovers et al. 1998): 



(B2) 



where the variable q here represents the column 
density of the electrons integrated along the ray 
path: 

9 = ^1 ^edl (B3) 

where ro is the classical electron radius and c 
is the speed of light. The negative phase delay, 
or phase advance given above is characteristic of 
propagation in a plasma. More appropriate for 
consideration here is the group delay, obtained 
by differentiating cj) = vAp^ with respect to fre- 
quency 

Agd = q/i^^ (B4) 

which constitutes a net additive delay that in- 
creases at lower frequencies. Numerically, this 
delay has the value 



1.34 X 10"^ Xe/i^^ sec 



(B5) 



where Xe is the electron column density in e 

_2 

cm 

Integrating the electron densities in the halo of 
a 10^*^ eV shower along a radial path from outside 
the shower down to a radius of 1 m and then back 
out again (to approximate the radar ray path) 
yields Xe — 4 x 10^ e~ cm~^, and the correspond- 
ing delay at 30 MHz is 0.6 ns, or about 6° of RF 
phase. A group delay of this order will have a 
negligible effect on the received power compared 
to the geometric phase factors. At lower frequen- 
cies, such at 10 MHz, the group delay is an order 
of magnitude larger, but the induced RF phase 
lag is still only 20°. Thus it appears that the 
refractive effects of the radial plasma halo can 
be ignored to first order for the EAS ionization 
columns considered here. 

Similar analysis of ray-path deviation along 
chords in the outer halo (due to the gradient in 
the plasma density) shows that the refraction is 
not significant until the chords approach within 
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several m of the EAS core for the highest en- 
ergy showers. Thus although it is important to 
include such effects in any final calibration of pa- 
rameters estimated by radar measurements, they 
do not appear to bear significantly on the de- 
tectability of the showers. 
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